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Chapitre 1 

Représentation des nombres en 
machine 

1.1 Arithmétique des calculateurs et Sources 
d’erreurs 

Si sophistiqué qu’il soit , un calculateur ne peut fournir que des réponses 
approximatives. Les approximations utilisées dépendent à la fois des con- 
traintes physiques (espace mémoire, vitesse de l’horloge...) et du choix des 
méthodes retenues par le concepteur du programme . (pour plus de détails 
sur le fonctionnement d’un ordinateur et la terminologie de base voir par 
exemple la page web htttp://www.co mmentcamarche.com 

Le but de ce chapitre est de prendre connaissance de l’impact de ces 
contraintes et de ces choix méthologiques. Dans certains cas il doit être pris 
en compte dans l’analyse des résultats dont une utilisation erronée pourrait 
être coûteuse. 

La première contrainte est que le système numérique de l’ordinateur est 
discret, c’est à dire qu’il ne comporte qu’un nombre fini de nombres; Il en 
découle que tous les calculs sont entachés d’erreurs. 

1.1.1 Evaluation de l’erreur 

Rappelons d’abord quelques notion de base ; 

Si X est une quantité à calculer et A"* la valeur calculée, on dit que : 

1. X — X* est l’erreur et \ E \—\ X — X* I est l’erreur absolue. 
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Exemple : 

Si A = 2.224 et A"* = 2.223 alors l’erreur absolue 

\ E \ = \ X — X* I = 2.224 - 2.223 = 0.001 


2 . Er = 


x-x* 


X r 


est l’erreur relative, X r 0. X r est une valeur de référence 
pour X. En général ,on prend X r = X. 

Exemple : 

Si X = 2.224 et A"* = 2.223 alors , si on prend X r = X, l’erreur relative 

0.001 


E r = 


* 

1 


X 

- X* I 

X r 


A 



2.224 


4.496 x 10 


-4 


Cependant, si X est la valeur d’une fonction F(t) avec a < t < b. on 
choisira parfois une valeur de référence globale pour toutes les valeurs 
de t. 

Exemple : 

SI X = sin(t) avec 0 < t < |, on pourra prendre 


A 


V2 

2 


sup sin(f). 

o<t<f 


En général , on ne connait pas le signe de l’erreur de sorte que l’on 
considère les erreurs absolues et les erreurs relatives absolues. 

Les opérations élémentaires propagent des erreurs. 

Dans la pratique, on considère que : 

1) L’erreur absolue sur une somme est la somme des erreurs absolues. 

2) L’erreur relative sur un produit ou un quotient est la somme des 
ereurs relatives. 

On peut estimer l’effet d’une erreur E sur l’argument x d’une fonction 
f(x ) au moyen de la dérivée de f(x). En effet f(x + E) ~ f(x) + Ef'(x) 

Exemple : 

Calculer la valeur de (llllllll) 2 

La valeur fournie par une petite calculatrice à cinq chiffres est 1, 2345xl0 14 
Mais la réponse exacte est 123456787654321. 
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La machine a donc tronqué le résultat à 5 chiffres et l’erreur absolue est 
de 6 * 19 9 . 

L’erreur relative est de 0.0005% . 

Cet exemple montre qu’il faut établir clairement l’objectif visé. 

Cet objectif est double ; 

1) Nous voulons un bon ordre de grandeur (ici 10 14 ) et avoir le maximum 
de décimales exactes, 

2) Ce maximum ne peut excéder la longueur des mots permis par la 
machine et dépend donc de la machine 


1.1.2 La mémoire de l’ordinateur : le stockage des 

nombres 

La mémoire d’un ordinateur est formée d’un certain nombre d’unités adess- 
ables appelées OCTETS . Un ordinateur moderne contient des millions voir 
des milliards d’octets. Les nombres sont stockés dans un ordinateur comme 
ENTIERS ou REELS. 

Les nombres entiers : 

Les nombres entiers sont ceux que l’on utilise d’habitude sauf que le plus 
grand nombre représentable dépend du nombre d’octets utilisés: 

-avec deux (2) octets, on peut représenter les entiers compris entre 

-32768 et 32767 

-avec quatre (4) octets on peut représenter!- les entiers compris entre 
-2147483648 et 2147483647 


Les nombres réels 

Dans la mémoire d’un ordinateur, les nombres réels sont représentés en no- 
tation flottante. 

Cette notation a été introduite pour garder une erreur relative à peu prés 
constante; quelque soit l’ordre de gandeur du nombre qu’on manipule. 

En notation flottante, un nombre a la forme: 
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x = ±Y x b' 


b est la base du système numérique utilisé 

Y est la mantisse : une suite de s entier yiy 2 ---y s avec y\ ^ 0 si x ^ 0 et 
0 < yi < (b - 1) 

e est l’exposant (un nombre entier relatif) 

La norme choisie est celle où la mantisse est comprise entre 0 et 1 et où 
le premier chiffre après la virgule est différent de zéro. 

Calcul de l’erreur 

Nous terminons ce chapitre en définissant les notions de troncature et 
d’arrondie. 

Exemple : 

En base 10, x — 1/15 = 0.066666666 

Dans le cas d’une représentation tronquée nous aurons , pour s = 5, 
fl(x) = 0.66666 *10" 1 . 

Remarquez comment nous avons modifié l’exposant afin de respecter la 
règle qui veut que le premier chiffre de la mantisse ne soit pas nul . 

Dans ce cas, l’erreur absolue A" — fl(X ) est de 6 x KL 7 '. L’erreur relative 
est de l’ordre de 10 -5 

Dans une représentation tronquée à s chiffres, l’erreur relative maximale 
est de l’ordre de 10 -s 

Dans une représentation arrondie, lorsque la première décimale négligée 
est supérieure à 5, on ajoute 1 à la dernière décimale conservée. 

Exemple : 

x = 1/15 = 0.066666666. 

Nous écrirons fl(x ) = 0.66667 x 10 -1 

L’erreur absolue serait alors 3.333 x 10 ~ 7 et l’erreur relative serait 5 x 
10“ 6 

En général, l’erreur relative dans une représentation arrondie à s chiffres 
est de 5 x 10' soit la moitié de celle d’une représentation tronquée. 
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1.2 Les régies de base du modèle 

Pour effectuer une opération sur deux nombres réels, on effectue l’opération 
sur leurs représentations flottantes et on prend ensuite la représentation flot- 
tante du résultat, 
l’addition flottante 


x®y = fl(fl(x) + fl(y)) 
la soustraction flottante 

xOy = fl((x) - fl(y )) 
la multiplication flottante 

x®y = fl{fl{x) x fl(y )) 

la division flottante 

x + y = fl(fl(x)/fl(y)) 


Chaque opération intermédiaire dans un calcul introduit une nouvelle 
erreur d’arrrondi ou de troncature. 

Dans la pratique, il faudra se souvenir du fait que deux expressions 
algébriquement équivalentes peuvent fournir des résultats différents et que 
l’ordre des opérations peut changer les résultats. 

Pour l’addition et la soustraction on ne peut effectuer ces 2 opérations 
que si les exposants sont les mêmes. On transforme le plus petit exposant et 
donc on ne respecte plus la régie voulant que le premier chiffre de la mantisse 
ne soit pas nul. 

Quelques remarques sur ce modèle: 

On constate une déviation importante par rapport aux lois habituelles de 
l’arithmétique. 

x + (y + z) peut être différent de (x + y) + z. 

Exemple : 

Pour 4 chiffres significatifs (s = A) on a : 

(. 1 + 0.0005) + 0.0005 = 1.000 
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et 


0.1 x 10 1 +0.5. x 10~ 3 = 0.1. x 10 1 +0. 00005. x 10 l = 
0.1 x 10 l +0.0000. x 10 1 = 0.1 x 10 1 


1 + (0.0005 + 0.0005) = 1.001 

Ainsi, l’addition flottante n’est, pas associative .(TD: Sommation d’une 
série à termes positifs) 

On constate aussi que si y est très petit par rapport à x, l’addition de x 
et y donnera seulement x. 

Exemple : 

L ’équation l + x = xax = 0 comme unique solution. Mais dans un 
système à 10 chiffres significatifs, elle aura une infinité de solutions (il suffit, 
de prendre | a: | — < 5 x 10~ n ) 

La distributivité de la multiplication par rapport à l’addition. 
Exemple : 

Considérons l’opération 

122 x ( 333 + 695) = ( 122 x 333 ) + (122 x 695) = 125416 

Si nous effectuons ces deux calculs en arithmétique à 3 chiffres (s = 3) 
et. arrondi, nous obtenons: 

122 x (333 + 695) = fl(122 ) x fl(1028) 

= 122 x 103 x 10 1 = fl(125660 ) = 126 x 10 3 
(122 x 333 ) + (122 x 695) = fl(40626) + fi(84190) 

406 x 10 2 +848 x 10 2 = fl(4 06 + 848) x 10 2 = fl(1254 x 10 2 ) = 125 x 10 3 

Donc la distributivité de la multiplication par rapport, à l’addition n’est 
pas respectée en arithmétique flottante. 



1.3 Propagation des erreurs. 


Une étude de la propagation des erreurs d’arrondi permattra d’expliquer ce 
phénomène. 

Soit à calculer e x à l’aide de son développement en série qui est convergent 
pour tout x : 


x x 

= 1+ IÏ + 2! + '" 


Il est évident que dans la pratique il est impossible d’effectuer la somma- 

/c 

tion d’une infinité de termes. On arrêtera donc lorsque le terme général U 
devient inférieur à 10 _ *(on a t digits). Pour x négatif on sait que le reste de 
la sérié est inférieur au premier terme négligé donc à 1Ü _/ (puisque la sérié 
est altérnée). 

Les calculs suivant sont fait sur ordinateur pour t = 14. 


X 

e x 

S 

1 

H 

O 

4.54. 10 -5 

4.54. 10 -5 

-15 

3.06. 10 -7 

3.05. 10 -7 

-20 

2.06. 10 -9 

— 1.55.10 -7 

-25 

1.39.10 -11 

1.87.10 -5 

-30 

9.36. 10 -14 

6.25. 10 -4 


On voit que pour x < 20 les résultats obtenus sont dépourvus de sens. 
L’explication de ce phénomène est la suivante: pour x = —30 les termes 
de la sérié vont en croissant jusqu’à (U- = 8.10 11 puis ils décroissent et 


æ 107 ~ 

107 ! 


9.19.10 -15 . 


L’erreur absolue sur le terme maximal est de (8.10 1 1 .HU 15 = 8.10 4 . Ainsi 
le résultat obtenu pour S représente uniquement l’accumulation des erreurs 
d’arrondi sur les termes de plus grand module de développement en sérié. 


1.3.1 Conditionnement et stabilité numérique. 

Le fait que certains nombres ne soient pas représentés de façon exacte dans 
un ordinateur entraine que l’introduction même de donnée d’un problème en 
machine modifie quelque peu le problème initial; Il se peut que cette petite 
variation des données entraine une variation importante des résultats. C’est 
la notion de conditionnement d’un problème. 

On dit qu’un problème est bien (ou mal) conditionné, si une petite varia- 
tion des données entraine une petite (une grande) variation sur les résultats. 
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Cette notion de conditionnement est liée au problème mathématique lui 
même et est indépendante de la méthode utilisée pour le résoudre. 

Une autre notion importante en pratique est celle de stabilité numérique. 
Un problème peut être bien conditionné et la méthode utilisée pour le ré- 
soudre peut être sujette à une propagation importante des erreurs numériques. 

Ces notions de conditionnement d’un problème et de stabilité numérique 
d’une méthode de résolution sont fondamentales en analyse numérique. Si un 
problème est mal conditionné alors la solution exacte du problèmetronqué ou 
arrondi à t digits pourra être très différente de la solution exacte du problème 
initial. Aucune méthode ne pourra rien; il faudra essayer de donner une autre 

formulation au problème. 

1 


X S. El Bernoussi, S. El Hajji et A. Sayah 
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Chapitre 2 

Résolution de f(x)=0 


2.1 Introduction 

Soit f une fonction numérique d’une variable réelle. 

On cherche les racines simples de l’équation 

(1) f(x) = 0 

La première étape consiste à isoler les racines, c’est à dire trouver un in- 
tervalle [a, b] dans lequel a est l’unique racine réelle de (1). On supposera que 
/ est continue et dérivable autant de fois que nécessaire dans cet intervalle. 

Pour trouver cet intervalle on aura besoin de quelques calculs prélim- 
inaires en utilisant soit le graphe des fonctions, soit (si la fonction / est 
continue dans [a, 6]) le théorème des valeurs intermédiaires en calculant /(a) 
et /(&) 

Si /(a) * f (b) < 0 / admet un nombre impair de racines dans 

[a, b] 

Si f(a)*f(b)> 0 / admet un nombre pair de racines 

Exemple : 

Soit, la fonction du grahe suivant : 

La fonction n’est pas définie pour x = ln( 2) et on a f'(x) — 1 + ^x1 2 ) 2 
donc f'(x) > 0 pour tout x. 

L’équation a donc 2 racines simples situées de chaque côté de ln( 2). 

On vérifie sans problème qu’une première racine appartient à [—1,0] et 
la deuxième à [1,2] 
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On supposera donc désormais avoir trouvé un intervalle [a, b] où / admet 
une unique racine simple et on supposera que / est définie, continue, et 
autant de fois continûment dérivable que nécessaire. 

Nous allons à présent définir la notion d’algorithme. 

Définition : Nous appellerons algoritnme toute méthode de résolution 
d’un problème donné. 

Pour tout problème, nous avons des données et des résultats. Les données 
sont appelées paramètres d’entrée (input) et les résultats paramètres de sortie 
(output). Ils constituent l’interface de l’algorithme (ou encore la partie visible 
de l’algorithme). 

Dans ce chapitre, nous désignerons par {p n } une suite de nombres réels . 

Il y a plusieurs façons de générer les termes d’une suite. En analyse 
numérique, on construit les suites à l’aide d’un procédé itératif appelé algo- 
rithme. 

Les algorithmes classiques que nous allons étudier sont les suivants: 

i) Méthode de la bissection 

ii) Méthode de Newton-Raphson 

iii) Méthode de la sécante 

iv) Méthode du point fixe. 

Le but de ce chapitre est de trouver des approximations de la solution de 
l’équation (1) avec une précision donnée et un nombre d’itérations maximum. 
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Afin de comparer ces différentes méthodes, nous allons introduire la no- 
tion d’ordre de convergence. 


2.2 Méthode de la bissection. 

Considérons une fonction f(x) quelconque, continue et cherchons p tel que 

f(p) = o. 

Nous supposons qu’on a localisé par tâtonnement un intervalle [a, 6] dans 
lequel la fonction change de signe (c.à.d. /(a) * f(b ) -< 0) on pose c = g ^, si 
/(a) * / (c) -< 0 on remplace b par c sinon on remplace a par c, et on continue 
cette operation jusqu’à ce qu’on trouve p avec la précision demandée. 

Algorithme de bissection (ou de dichotomie) 

But : Donner une fonction continue f(x) et un intervalle [a, b] pour 
lequel /(a) et f{b) sont de signes contraires, trouver une approximation 
de la solution de f(x) = 0 dans cet intervalle; en construisant une suite 
d’intervalles ([a n ,b n }) n contenant cette racine et têts que a n ou b n est le mi- 
lieu de l’intervalle [a n _i, b n -i\. 

Entrées : a, b les extrémités de l'intervalle 
e la précision désirée 
N 0 le nombre maximal d’itérations 
Sortie : la valeur approchée de la solution de f(p) = 0 

Etape 0: Si /(a) = 0 imprimer la solution est a, aller à l’étape 9 
Si f{b) =0 imprimer la solution est 6, aller à l’étape 9 

Etape 1: 

si /(&) * /(a) > 0 

alors imprimer (il n’y a pas de changement de signe) 
aller à l’étape 9 
Etape 2: poser N = 1 

Etape 3: 

Tant que N < N 0 , faire les étapes 4 à 7 
Etape4: poser p = 

Etape 5:Si f(p ) = 0 ou ^ < e 
Alors imprimer p 
aller à l’étape 9 
Etape 6: poser N = N + 1 
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Etape 7: Si /(a) * f(p) > 0 
alors poser a = p 
sinon poser b = p 

Etape 8: Imprimer après N 0 itérations l’approximation obtenue est p et 
l’erreur maximale est 

Etape 9: Fin 


2.3 Méthode de Newton-Raphson: 

Le principe consiste à construire une suite (x n ) m telle que x n+ \ soit 
l’intersection de la tangente à la courbe de / au point (x n , f(x n )) avec l’axe 
horizontal. 



Figure 2.2: Méthode de Newton pour f(x) = log(x),x 0 = 2. 


On a: 


donc 


M = Oo, f(x o)), B = (x i, 0) G axe(Ox) 
[ A et B G D : y = ax + b 

f f(x 0 ) = ax 0 + b \ a = f'(x 0 ) 

\ 0 = ax 1 + b y x 1 = x 0 - jj&L 


Algorithme de Newton-Raphson. 
But: Trouver une solution de f(x) = 0 
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Entrées: une approximation initiale po 
£ (la précision désirée) 

Nq (le nombre maximum d’itérations) 

Sortie: valeur approchée de p ou un message d’échec 

Etapel : N — 1 

Etape 2: Tant que N < N 0 , faire les étapes 3 à 6. 

Etape 3: Poser p = p 0 — 

Etape 4: Si | p — po |< e alors imprimer p 
aller à l’étape 8. 

Etape 5: Poser N — N + 1. 

Etape 6: Poser p 0 = p. 

Etape 7: Imprimer la méthode a échoué après N itérations. 

Etape 8: Fin. 


2.4 Méthode de la sécante 

La méthode de Newton- Raphson suppose le calcul de f'(p ) à chaque étape. 
Il se peut qu’on ne dispose pas d’un programme permettant de calculer sys- 
tématiquement f . 

L’algorithme suivant peut être considéré comme une approximation de la 
méthode de Newton. 

Le principe consiste à construire une suite (x n ) n à l’aide de la formule 
obtenue en remplaçant dans la méthode de Newton f'(p n ) par ^ P, ^Z • 
Ainsi au lieu d’utiliser la tangente au point p n nous allons utiliser la sécante 
passant par les points d’abscisses p n et p n -i pour en déduire p n+ \. Ce dernier 
est obtenu comme intersection de la sécante passant par les points d’abscisse 
Pn et Pn- 1 et de l’axe des abscisses. 

L’équation de la sécante s’écrit : 
s(x) = f(p n ) + (x- p n ) 

Si s(p n+ i) = 0 , on en déduit: 

Pn+l — Pn / (Pn) /(pC)-/(p,!- 1) 

Algorithme de la sécante: 

But:Trouver une solution de f(x) = 0 
Entrées: deux approximations initiales po et [>\ 
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Figure 2.3: f(x) = x 3 — 1 


£ (la précision désirée) 

N 0 (le nombre maximum d’itérations) 

Sortieda valeur approchée de p ou uu message d’échec 
Etape 1: poser N — 1 
(lo = fipo) 

Qi = f(Pi) 

Etape 2: Tant que N < N 0 + 1, faire les étapes 3 à 6 
Etape 3: poser p = pi — 

Etape 4: Si | p — p± \ < e alors imprimer p 
aller à l’étape 8 
Etape 5: Poser N = N + 1 
Etape 6: Poser p 0 = p 1 
% = Qi 
Pi = P 
Qi = f(P ) 

Etape 7: Imprimer la méthode a échoué après Nq itérations 

Etape 8: Fin. 


2.5 Méthode du point fixe 

Nous pouvons observer que la méthode de Newton peut s’interpréter comme 
Pn+ 1 = g(Pn) OÙ 
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g(x) = x — ( jt^T) )• Maintenant , si la fonction g(x) est continue et si 
l’algorithme converge (c.à.d. p n — > p), on tire de p n+ i = g(p n ) que p satisfait 
l’équation p = g(p) ; on dit que p est un point fixe de g. 

On peut toujours transformer un problème du type f(x) = 0 en un prob- 
lème de la forme x = g(x) et ce d’une infinité de façons. 

Par exemple 

x 2 - 2 = 0 

ou x = 2/x 

ou x = x 2 + x — 2 

ou x = a(x 2 — 2) + x 

Il faut toutefois noter que ce type de transformations introduisent des 
solutions ’parasites’. 

Par exemple : résoudre l/x — a ou encore x = 2x — ax 2 

On voit que 0 est racine de la deuxième équation mais pas de la première. 

Algorithme du point fixe 
But: trouver une solution de g{x) = x 
Entrées: une approximation initiale /y 
c(la précision désirée) 

Nq le nombre maximale d’itérations 
Sortie: valeur approchée de p ou un message d’échec 
Etape 1: poser N = 1 

Etape 2: Tant que N < No, faire les étapes 3 à 6 
Etape 3: poser p = g(po) 

Etape 4: Si j p — p 0 |< e 
alors imprimer p 
aller à l’étape 8 
Etape 5: poser n — n + 1 
Etape 6: poser p 0 = p 

Etape 7: Imprimer (la méthode a échoué aprrès N 0 itérations) 

Etape 8 : Fin. 
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2.6 Convergence et ordre de convergence. 

Définition: Soit D une partie de R. et F une application de D dans D. On 
dit que la fonction F est contractante si 

Vx, y G D , 3k G [0, 1[ tel que 

| F(x ) — F (y) | < k \ x — y \ . 

k est le coéfficient de contraction ou de Lipschitz de F. 

Théorème: Considérons le segment S = [po — a,po + a] C D] si F est 
contractante sur S et si | F(p 0 ) — p 0 |< (1 — k) a, alors l’itération p n+ \ = 
F{p n ) de point initial p 0 , converge vers l’unique point fixe p G S de F. 

Théorème: Convergence locale. 

Si F est différentiable au voisinage d’un point fixe p et si | F\p) |< 1 
alors : 

3V voisinage de p tels que po G V et p n +i = F{jp n ) converge vers p. 
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2.6.1 Interprétation graphique. 



Figure 2.4: F{x) = — x 3 + |x; |F'(o;)| < 1, convergence. 


On voit graphiquement que | F'(p) |< 1, et par conséquent les itérations 
convergent vers le point fixe, p est un point fixe attractif. Par contre si 
| F' (p) |> 1 pas de convergence vers le point fixe, p est un point fixe répulsif. 



F(x) = X* — x, |F'(16)| > 1, l’itération diverge. 

Remarque: Un point fixe répulsif pour une méthode devient attractif 
pour une autre. 
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2.6.2 Ordre de convergence. 

La convergence de l’itération p n +i = F(p n ) vers le point fixe peut se 
faire plus ou moins vite. 

Définition : Considérons une suite {p n } convergeant vers p et posons 
e n = Pn - P- 

On dit dans le cas où j j converge, que la suite p n converge linéaire- 
ment vers p ou encore que la méthode est du premier ordre. 

Si on a j (e £n ^ k j converge, alors la convergence est dite d’ordre k. 

Exemple : 

La méthode de Newton pour résoudre l’équation f(x) = 0 est une méth- 
ode de type point fixe avec F(x) = x — jrp^- Si x* est racine simple de 
f(x) = 0, alors f'(x*) ^ 0 et il existe un voisinage V de x* tel que pour tout 
Po G V, la suite (p n ) n converge vers x* et l’ordre de convergence est 2. 

(en effet F'(x) = 1 — ^ ^ ^ F'(x*) = 0. Ainsi d’après 

le théorème précédent la méthode de Newton converge. Pour déterminer 
l’ordre de convergence on utilise la formule de Tylors en x* : F(x) = F(x*) + 
F'(x*)(x - x*) + F"(9 x) { ^=££). 

2.7 Exercices 
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Série f(x) = 0 


Exercices 1 

Résoudre à l’aide de la méthode de bisection tanx — x = 0 dans l’intervalle 


[4; 4.7], 


Exercice 2 

On considère l’équation 

(1) e x - Ax = 0 

1) Déterminer le nombre et la position approximative des racines de (1) 
situées dans .x > 0 

2) Utiliser l’algorithme de dissection pour déterminer la plus petite de 
ces racines à £ près, (par exemple HR 7 ) 

3) Sans faire d’itérations, déterminer combien vous devriez en faire pour 
calculer la plus grande racine à l’aide de la bissection avec une précision de 
1CR 8 , si l’intervalle de départ est [2; 2, 5] 

Exercice 3 

Ecrire un algorithme pour calculer par la méthode de Newton la racine 
K-ième d’un nombre. 

Quelle est la valeur de s = \J 2 + y^2 + \/2 + ? 

Suggestion: écrire .p n+ i=G(p n ),p o =0 Quel est l’ordre de convergence ? 

Exercice 4 

Ecrire 3 méthodes itératives pour la résolution de a: 3 — x — 1 =0 et vérifier 
expérimentalement leur convergence avec xo = 1,5. Trouver à HR 6 près la 
racine comprise entre 1 et 2. Connaissant la valeur de cette racine, calculer 
l'ordre de convergence de vos 3 méthodes. Ce résultat coincide-t-il avec 
l'expérience? 

Exercice 5 

Résoudre x 2 -l=0 en utilisant la méthode de la sécante avec x 0 = —3 et 
Xi = 5/3. Qu’arrivera-t-il si on choisit et x 0 = 5/3 et X\ = —3? Expliquez. 

î 


3 S. El Bernoussi, S. El Hajji et A. Sayah 
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Chapitre 3 
Algèbre linéaire 


3.1 Introduction 

Un sytème linéaire s’écrit sous la forme : 

(1) Ax = b 

où A est une matrice rixn à coefficients réels, b G M™ et x G M n . 

La résolution de grands systèmes linéaires (et non linéaires) est pratique 
courante de nos jours. Elle apparait dans tous les domaines où l’on s’intéresse 
à la résolution numérique d’équations aux dérivées partielles. 

Il existe plusieurs packages (linpack, eispack, ..), logiciels (Maple et Mat- 
lab) et programmes (http: / / www.netlib.com, numerical recipes, NAG, IMSL, 
...) de base pour le résoudre. 

Le choix de la méthode dépend fortement du type (forme) de la matrice. 
Les méthodes de résolution sont de deux types : 

Les méthodes directes : Une méthode est dite directe si elle permet d’obtenir 
la solution en un nombre fini d’opérations. 

Les méthodes itératives : Une méthode est dite itérative si elle permet de 
construire une suite (. x n ) n qui converge vers la solution. 

Dans ce chapitre nous allons : 

1. Rapeler des notions et notations de base relatives aux systèmes linéaires 
et aux matrices 
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2. Etudier une méthode directe : la méthode de Gauss. 

3. Etudier la décomposition (factorisation) LU. 

4. Etudier des applications : Inverse de matrices,... 


3.2 Rappels sur les systèmes linéaires 

Un système de n équations linéaires à n inconnues peut toujours s’écrire sous 
la forme : 

(1) Ax = b 

où A est une matrice (a^) et x et b sont des vecteurs colonnes de dimension 
n. 

Si la matrice A est inversible alors le système linéaire (1) admet une 
unique solution x = /I 1 b où A' 1 est la matrice inverse de A. 

Ainsi théoriquement le problème revient à calculer A -1 ? Mais en pratique 
ce calcul est difficile. 

Il existe plusieurs méthodes classiques pour résoudre (1) sans calculer 

A- 1 . 

Pour cela on va considérer le cas simple suivant : 


( 1 ) 


f x + 2y = 5 
\ 2x + y = 4 


i) La méthode de Cramer consiste à calculer la solution en calculant des 
déterminants. 



5 2 


1 5 


On a: x — ■ 

4 1 

- = zf = 1 et y = 

2 4 

. — Ai — 2 

1 2 

1 2 


2 1 


2 1 


ii) La méthode de substitution (ou d 

/élimination) consiste 


le système (1). 


( 1 ) 


x + 2y = 5 
2x + y = 4 


x = —2 y + 5 
2x + y = 4 


x = —2 y + 5 
2(-2y + 5) + y = 4 
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J x = —2 y + 5 J x = —2 y + 5 j x — 1 

^ l 32/ = 6 ^ \ y = 2 ^ \ y = 2 

Peut-on généraliser ces méthodes pour un système de n équations avec 

n G N? 

Théoriquement OUI mais en pratique cela va nécessiter beaucoup de cal- 
culs et de techniques. 


3.3 Méthode Gauss 


La méthode de résolution la plus étudiée (et une des plus employées) s’appelle 
méthode d’élimination de Gauss. 

L’idée de base de cette méthode consiste à transformer le système linéaire 
(1) en un problème que l’on sait résoudre. 

Si la matrice A = D avec D une matrice diagonale, alors on sait résoudre 

( 1 )- 

Mais toute matrice n’est pas diagonalisable. 

Si la matrice A = U (ou L ) avec U (ou T) une matrice triangulaire 
supérieure ( ou inférieure) alors on sait résoudre (1). 

Problème : Comment tranformer une matrice en une matrice triangulaire 
inférieure ou supérieure ? 

La méthode de substitution (d’élimination) répond à cette question mais 
elle n’est pas automatique. 

La méthode d’élimination de Gauss a pour but de remplacer le sys- 
tème (1) par un système triangulaire possédant la même solution. Son 
principe s’apparente à celui de la méthode de substitution (d’élimination) 
mais (comme on le verra ci dessous), il est plus simple à automatiser. 

Regardons son fonctionnement sur l’exemple suivant cas n — 3: 

On pose A = (aÿ)ij=i,3 X = (xj^i^et b = (&i)i= lj3 de telle sorte que 
AX = b s’écrit sous la forme : 

{ aiixi + <112X2 + 0*13X3 = b\ 
a, 21X1 + 0,22X2 + 023X3 = 62 
031X1 + 032X2 + 033X3 = 63 


ou encore sous la forme dite augmentée 


(A b) 


0 11 Oi2 0\3 b 1 

021 °22 °23 ^2 

®31 a 32 Û33 63 
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On suppose que an ^ 0 , par élimination, on obtient : 


{A 1 h) 


«Il «12 «13 W \ 

0 a 22 a 2 3 b 2 J 

0 « 3 2 «33 b 3 J 


On va illustrer la méthode de Gauss sans passer par le sytéme augmenter 


On a : 


{ «11«U + «12*^2 + «13^3 = b\ (/l) 
«21«h + d 22 X 2 + «23«'3 = &2 (G) 

«31«h + «32^2 + «33^3 = &3 (h) 

On note par (U) la i eme équation du système précèdent. 
On suppose que an ^ 0 , 

On pose : 

(l 2 ) = aufy — a2i(/ / 1 ) 
et 

(I3) — « 11 (^ 3 ) — « 31 (^ 1 ) 

Alors (1) s’écrit 


( 2 ) 


aiPi + a 12 x 2 + 013X3 = 61 (11) 

a 22 x 2 + a' 23 x 3 = b' 2 (V 2 ) 

o! 32 x 2 + a' 3 3X3 = b' 3 (V 3 ) 


On suppose que a 22 ^ 0, 

On pose : 

(l 3 ) = a 22 (l 3) — «32(^2) 

Alors (2) s’écrit 

{ «11X1 + 032X2 + 013X3 = b \ ( 11 ) 

a! 22 x 2 + 033X3 = U 2 (l' 2 ) 

«33^3 = &3 (O 

Remarque : 

i) Les termes diagonaux à chaque étape sont appelés les pivots, 
ii) Si un pivot a ü est nul on change de ligne (on permute) de i à n (pivotage 
partiel) 
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iii) Cette méthode se généralise assez facilement bien qu’il faut être pru- 
dent avec le choix du pivot. En pratique, il faut éviter de prendre des pivots 
"trop" petits. 

Exemple : Sur l’importance du pivot 
1) On considère le système : 


( 1 ) 


f x + y = 2 
\ 10 _4 a: + y = 1 


Calculer la solution exacte de ce système. 

2) Calculer la solution pour s = 3 avec troncature des systèmes 


( 1 ) 


f x + y = 2 
\ 10” 4 a; + y = 1 


et 


f 10” 4 a; + y = 1 
[ x + y — 2 


Remarque : L’algorithme de Gauss est une méthode systématique de résolu- 
tion de systèmes d’équations comportant un nombre quelconque d’inconnues. 

Dans le cas où tous les pivots sont non nuis i.e. an ^ 0, l’algorithme: 
Elimination de Gauss s’écrit : 

Partie 1: Réduction à la forme triangulaire (ou élimination de Gauss) 
Entrée A et b 

Sortie A = U (forme triangulaire), et b. 

Pour j = 1, ..., (n — 1) 

Pour i — j + 1, ..., n 

7 , a ij 

".il 

Pour k — j + 1, ..., n 
& ik * ^ik lijüjk 

Fin 

bj — hjbj 

Fin 

Fin 

Sortie A = U (forme triangulaire), et b 

Cette partie s’écrit sous la forme (algorithmique) 

Etape 1 : Poser j = 1 

Etape 2: Tant que j < n — 1 faire 

Etape 3: Si a.jj = 0 afficher ’pivot nul’ aller à étape 14, sinon 


26 



Étape 4: Poser i — j + 1 
Étape 5: Tant que i < n faire 
Étape 6: Uj = < ^ L 

Étape 7: Si = 0; aller à l’étape 12. 

Étape 8: Poser k — j + 1 
Étape 9: Tant que k < n faire 

Étape 10: — hjdjk, k = k + 1; aller à l’étape 9. 

Étape 11: bi — bi — hjbj ; 

Étape 12: poser i — i + 1; Aller à l’étape 5. 

Étape 13: j — j + 1; Aller à l’étape 2. 

Étape 14: Fin. 

Remarques:Les éléments sous la diagonale principale de la nouvelle matrice 
obtenue sont nuis. Comme ils n’interviennent pas dans la résolution du 
système triangulaire formé, il est inutile que l’algorithme leur assigne cette 
valeur nulle. 

Exemple : On considère le système linéaire : 

{ x + y + 3f = 4 
2x + y — z + t = 1 
3x — y — z + 2t = — 3 
— x + 2y + 3z — t = 4 

qui s’écrit encore: 


/ 1 1 0 3 \ 


( x \ 


( 4 \ 

2 1-11 


y 


i 

3-1-1 2 


z 


-3 

\ -1 2 3 -1 ) 


V t y 1 


l 4 / 


Nous appliquons l’algorithme à 
augmentée. 

Nous obtenons 

A 

1 1 0 
2 1-1 
3 -1 -1 
-12 3 


exemple en travaillant sur la matrice 


b 

3 . 4 

1 . 1 
2 . -3 
-1 . 4 


notre 
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1 

1 

0 

3 . 

4 

0 -1 

-1 

-5 . 

-7 

1 

O 

-1 

-7 . 

-15 

0 

3 

3 

2 . 

8 

1 

1 

0 

3 . 

4 

O 

1 

-1 

-5 . 

-7 

0 

0 

3 

13 . 

13 

0 

0 

0 

-13 . 

-13 


Que l’on peut écrire sous la forme : 

{ x + y + 3t = 4 
—y — z — 5t = —7 
3z + 13/ = 13 ’ 

— 13t = -13 

Notons que l’étape j — 3 nous donnerait / 43 = 0. 

Nous avons maintenant un système triangulaire à résoudre. 

Partie 2 : Remontée triangulaire 

Entrée A, b avec A matrice triangulaire supérieure 
Sortie x solution du sytème Ax = b 

• Étape 1: x n = — 

• Étape 2: Pour i — n — 1, n — 2, 1 faire: 

1 n 

Xi ( bi ^ ^ e, /j .v j ) 

• ■ , -, 

3=i+l 

En appliquant cet algorithme à notre exemple, nous obtenons x = (—1, 2, 0, 1). 

Remarque: 

1. Dans la pratique le test (3) de l’algorithme d’élimination de Gauss ne 
conduit pas à l’arrêt. En fait, si le pivot est nul, on cherche, dans 
la même colonne, un élément d’indice plus grand non nul, puis on 
échange les lignes correspondantes. Si ceci est impossible, le système 
est singulier. 
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2. On est parfois amené, pour des raisons de stabilité numérique, à ef- 
fectuer des échanges de lignes même si le test (3) est négatif (c’est à 
dire que le pivot est non nul). Ceci conduit à des stratégies dites de 
pivot que nous n’étudierons pas ici. 

Exemple : Résolution du système suivant : 

( 2x + 6y + lOz = 0 / 2 6 10 \ / x \ / 0 \ 

l x + 3y + 3z = 2 ^ 1 3 3 y = 2 

{ 3x + Uy + 2Sz = -8 \ 3 14 28 / \ z / \ -8 J 

/ 2 6 10 0 \ / 2 6 10 0 \ / 2 6 10 0 \ 

1 3 3 2 J | 0 0 —4 4 U 0 5 13 -8 

\ 3 14 28 -8 J \ 0 5 13 -8 / \ 0 0 -4 4 J 

En utilisant la rementé on trouve: 

| y = i(-8 — 13 x (-1)) = 1 => x’ = ( 2 ) 

( x= |(— 6 x 1 - 10 x (-1)) =2 \ -1 / 


3. Méthode de Gauss avec normalisation :Elle consiste à normaliser le 
pivot: 

On a : 

{ auXi + 012X2 + O13X3 = 61 ( h ) 

O21X1 + 022X2 + a,23 x 3 = b-2 (I2) 

031X1 + 032X2 + 033X3 = 63 ( Iz ) 


On note par (Zj) la i eme équation du système précèdent. 
On suppose que an ^ 0, 

(/ 1 ) s’écrit : xi + ^x 2 + ^x 3 = £ (Zi) 

Si on pose : 

{l' 2 ) = (h) - a 21 (l[) 

et 

(Zi) = (Z 3 )-fl 31 (Zi) 

Alors (1) s’écrit 

[ *1 + S 1 * + Sî** = & «) 

(2) < a 22 x 2 + a 23 x 3 = ^2 ( k ) 

l a 32 X 2 + a 33 X 3 = (Z s ) 
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On suppose que a 22 0) 

(l' 2 ) s’écrit x 2 + 023X3 = b '2 (l") 

Si on pose : 

{%) = &)-<&(%) 

si «"g 7^ 0 on pose (l™) x 3 = ph 

Alors (2) s’écrit 

f + gp* + gp* = £ K) 

(2) ) X 2 + <4,13 = K ® 

1 13 = i <'") 

1. Cette statégie est très utile pour calculer l’inverse d’une matrice. 

Nous pouvons nous demander s’il existe une relation entre la matrice de 
départ et la matrice triangulaire obtenue. Ce lien existe. 


3.4 Factorisation LU 


Matriciellement la méthode de Gauss consiste à multiplier la matrice A par 
la matrice Li de telle sorte que l’on ait : 


A 1 = L X A 


Li = 


1 

«21 

«11 

«31 

«11 


0 0 
1 0 
0 1 


On suppose que a 22 7^ 0, donc on cherche L 2 de telle sorte que 


A 2 — L 2 A± — 


=>■ l 2 = 


(1 

0 

v° 


an 

® 1 2 

«13 

0 

a 22 

a 

/ 

23 

0 

0 

a 

// 

33 


0 0 

\ 



1 0 




a 32 ^ 




a 22 

/ 
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Ainsi on a : A 2 = L 2 A 1 — U et A 2 est une matrice triangulaire supérieure. 
De plus si on pose A 0 = A alors U = A 2 = L 2 L\A 0 c’est à dire que 

U = L 2 L x A ^A = L 2 l L 2 l U 


On a Li et L 2 sont des matrices inversibles et triangulaires inférieures 
donc L 2 * L | est une matrice inversible et triangulaire inférieure. 

De même L = L^ 1 * L- 2 1 est une matrice inversible et triangulaire in- 
férieure 

Donc A = LU. 

Ainsi le système linéaire (1) AX = b s’écrit 


LUX = bo 


Ly = b avec L matrice triangulaire inférieure 
U X = y avec U matrice triangulaire supérieure 


En conclusion ( à admettre) la méthode de Gauss revient à décomposer la 
matrice A en un produit de deux (2) matrices triangulaires l’une supérieure 
U et l’autre inférieure L. 

Avec : 


L 


( 1 

hi 1 

^31 ^32 1 0 


\ 


y Inl 


l 


n,n — 1 


1 / 


où lij est défini à l’étape (6) de l’algorithme d’élimination et (si l’algorithme 
d’élimination n’exige pas d’échange de lignes). 


Nous ne démontrerons pas cette proposition. Nous nous contenterons de la 

vérifier sur notre exemple. 

Exemple : 


( 

1 

1 

0 

3 ^ 


( 

1 

0 

0 

0 ^ 


( 1 

1 

0 

3 \ 


2 

1 

-1 

1 



2 

1 

0 

0 


0 

-1 

-1 

-5 


3 

-1 

-1 

2 



3 

4 

1 

0 


0 

0 

3 

13 

\ 

-1 

2 

3 



V 

-1 

-3 

0 



^0 

0 

0 

- 13 / 
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Il y a une classe importante de matrices pour lesquelles l’élimination 
peut toujours s’opérer sans échange de lignes (i.e. le pivot a,jj n’est jamais 
nul pendant l’algorithme d’élimination). Ce sont les matrices à diagonale 
strictement dominante. 

Définition: Une matrice A est dite à diagonale strictement dominante si 
pour tout i — 1, 2, ...., n , on a : 


> £ 


IJ I 


est vérifiée. 

Remarque : Si la matrice est à diagonale strictement dominante alors elle 
est inversible. 


3.4.1 Appplications de la Factorisation LU 

Si l’on doit résoudre souvent un système où seul le membre de droite change 
ou son veut calculer l’inverse d’une matrice, il y a intérêt à effectuer la ré- 
duction à la forme triangulaire une fois pour toutes. 

En effet, si A — LU on peut résoudre: Ax = b eu résolvant Lz = b et 
U x = z. On a : 


(1) Ax — b 


(2) Lz = b 

(3) Ux = z 


Dans ce cas Ax = LUx = L(Ux ) = Lz = b . 

Les systèmes (2) et (3) étant triangulaires, la résolution ne nécessite que 
l’exécution d’une remontée et d’une descente triangulaire. 


Exemple : 

( 1 

0 

0 

A = LU = 

2 

1 

0 

3 

4 

1 


W 

-3 

0 


/ 4 \ 


\ 

( 1 

1 

0 

3 \ 


0 

-1 

-1 

-5 


0 

0 

3 

13 

/ 

^0 

0 

0 

"13/ 


tème Ax = 


1 

-3 


, on résoud le sys- 


V 4 ) 
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Lz — b =>• 


Ux = Z =>• 


^i = 4 

2^i + £2 = 1 ^ _ 

Szi + 4^2 + ^3 = 3 

— Z\ — 3z 2 “I - Z 4 . — 4 
—13x4 = —13 
3 x 3 + 13x4 = 13 
— x 2 — x 3 — 5 x 4 = —7 ^ ' r 
Xi + x 2 + 3 x 4 = 4 



3.5 Mesure des erreurs 

L’utilisation d’un calculateur pour implanter les algorithmes étudiés conduira 
inévitablement à des erreurs. Pour mesurer celles-ci, nous devons mesurer la 
distance entre le vecteur représentant la solution exacte x = (xi, x n ) et le 
vecteur x = (âq, ..., x n ) représentant la solution approchée. Nous pouvons, 
pour ce faire, utiliser la "longueur" usuelle de R n i.e.: 

n 

Iklh = 

1 

pourtant, dans la pratique on lui préfère souvent la longueur 

llxlloo = max Ixjl 

1 <i<n 

Par exemple si x = (1, —7, 2, 4) alors HxHoo = 7 . 

Exemple : 

Si x = (1, 1, 1, 1) alors HxH^ = 1 si x = (1.01, 1.1, 1, 1), on a 

||x — xlloo = 0.1 

Considérons alors le système 

/ 10 7 8 7 \ / xi \ / 32 \ 

7 8 6 5 x 2 _ 23 

8 6 10 9 x 3 “ 33 

7 5 9 10 / \ x 4 / \ 31 y 

dont la solution exacte est x = (1, 1, 1, 1) . 
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Si dans le membre de droite nous remplaçons b par: 


b= (32.06; 22.87; 33.07; 30.89) 


lions obtenons 

x = (9.19; -12.59, 4.49, -1.09) 
C’est-à-dire qu’une erreur relative de l’ordre de: 

= 3* 10” 1 

sur b a entraîné une erreur relative de l’ordre de 

= 13.52 

|| "^ Il oo 



sur la solution. 

Nous devons donc soupçonner que l’application de l’arithmétique finie à 
la résolution d’un tel système serait désastreuse. L’étude de cette question 
dépasse le cadre de ce programme. 


3.6 Exercices 
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Série Ax = b 


Exercice I - 

1) On considère le système linéaire : 


( 1 ) 


( 1 6 

V 1.0001 5 )\y) V 6 - 2 * * 0005 


Déterminer la solution X de ce système. 

2) Dans le système précédent, on remplace 6.0005 par 6, déterminer la 
solution A"* de ce nouveau système notée (2). 

3) Calculer les erreurs relatives sur les données et sur les résultats. 

4) Conclusion. 


Exercice II - 

Résoudre le système linéaire (1): 


{ x + 2y + 3z = 1 
2x + 6y + 10 z = 0 
3a; + lAy + 28 z = —8 

1) Par Gauss Classique 

2) Par Gauss avec pivotage partiel 

3) Par Gauss avec pivotage et mise à l’échelle (i.e. a ü = 1). 

Exercice III - 

1) En arithmétique flottante avec 2 chiffres significatifs ( s = 2 (s est le 
nombre de digits)) et arrondi, résoudre par élimination de Gauss, les systèmes 
linéaires (1) et (2). 


( 1 ) 


0.0001a; + y = 3 
x + 2y = 5 


2) Conclusion 


( 2 ) 


x + 2y = 5 
0.0001a; + y = 3 


Exercice IV - 

Soit la matrice A 


30 -20 -10 \ / 1 \ 

-20 55 -10 et b = 5 

-10 -10 50 / \ 2 J 
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1) Ecrire la matrice A sous la forme LU i.e. trouver LetU (sans pivotage) 
avec L matrice triangulaire Inférieure et U triangulaire supérieur. 

2) En déduire le déterminant de A 

3) Résoudre par Factorisation LU. le système linéaire AX = b 


Exercice VI - 

/ 1 -1 2 
Soit la matrice A — I — 2 1 1 

\ -1 2 1 

1) Ecrire la matrice A sous la forme LU i.e. trouver L et U avec L matrice 
triangulaire Inférieure et U triangulaire supérieur. 

2) Utiliser 1) pour calculer le déterminant de la matrice A. 

3) Utiliser 1) pour calculer l’inverse de la matrice A. 


X S. El Bernoussi, S. El Hajji et A. Sayah 
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Chapitre 4 

Interpolation polynomiale 


4.1 Introduction 

Nous abordons dans ce chapitre un nouveau type de problème, faisant inter- 
venir la notion d’approximation d’une fonction. 

Cette notion a déjà été rencontrée dans les cours d’analyse. 

Exemples : 

1) D’après la Formule de Taylor à l’ordre 5 de la fonction sin(x), on a : 

3 5 6 

VT G Voisi 0), sin(x) ~ x — — + — + sin®(£) — - où £ G Voisin) 

3! 5! 6! 

On a tronqué la formule de Taylor après l’ordre N (ici 5), on obtient : 
au voisinage de 0, une approximation de sin(x) par un polynôme de 
degré N (ici 5). 

l’erreur commise serait de l'ordre de sin ( ' 6 )(£)|y- où £ G Voisin) 

Ainsi avec ce type d’approximation, on a : 

• Si N = 3, sin(O.l) = (0.1) - ^ = g. 983 3 x 10” 2 

• Si N = 5, sin(O.l) = 0.1 - ^ = 9.983 3 x 10” 2 

Avec le logiciel Maple on a : sin(O.l) = 9. 983 3 x 10“ 2 

2) Avec les cours d’analyse / et 17, on ne connait pas d’expression ex- 
plicite de / = Jq e~ x2 dx 

Cependant d’après : 
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• La formule du trapèze I = f^e x2 dx ~ IM±ÎP1 — 1+ ^ 1 = 0.683 94 

• La formule de Simpson : / = e~ x2 dx ~ | [/(O) + 4/(|) + /(l)] = 

|(1 + 4e" ï Le" 1 ) = 0.74718 


En utilisant la méthode des trapèzes et en subdivisant (partageant) le 


segment (intervalle) [0, 1] en 10 intervalles égaux, on a : I = f* ■ 
à^ + èEhe-^ + è = 0.74621 


J dx ~ 


En utilisant la méthode de Simpson et en subdivisant (partageant) le 


segment (intervalle) [0, 1] en 10 intervalles égaux, on a : / = fi 1 


" dx ~ 


à e 1 + è Sti e 3X5 ' 2 + â EL ex P (- ü* - è) 2 ) 


+ ^ = 0.746 82 


Avec le Logiciel Maple, on a : e x2 dx = \\f^erf (1) = 0.746 82 

NB : erf() est "The Error Function". Elle est définie pour tout x par : 
er f ( x ) = 77 fo e_<2rft - 

Donc l’erreur relative ( la qualité de l’approximation) dépend du type 
d’approximation choisie. 

On ne connait pas à ce niveau du cours l’expression explicite de l’erreur. 

La notion d’approximation d’une fonction consiste à remplacer un prob- 
lème donné par un problème voisin (un problème majeur en analyse numérique) 

La question fondamentale serais de savoir la qualité de cette approxima- 
tion i.e. la solution (du problème approché) obtenue est -elle aussi voisine 
que l’on veut de la solution du problème initial. 

Remarque : En pratique la fonction / est connue explicitement, ou seule- 
ment par ses valeurs en quelques points. 

La notion d’interpolation polynomiale est la façon la plus simple d’obtenir 
une telle approximation. 

Théorème : (à admettre) 

Soit / une fonction continue dans [a, b] C IR, alors pour tout e > 0 
donné, il existe un polynôme P n de degré n tel que 


Max | f(x) — P n (x) | < e 

x£[a,b] 
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Ce théorème ne permet pas de construire (de déterminer explicitement) 
le polynôme P n . Il existe cependant un certain nombre de techniques (algo- 
rithmes) qui le permettent : 

1. L’interpolation polynomiale: Elle est la plus classique et est un outil 
pour la construction des méthodes d’intégration numérique ou des méth- 
odes d’approximation des équations différentielles. 

Remarque : Pour les équations aux dérivées partielles, la méthode 
des éléments finis, un des outils de base de l’ingénierie moderne, utilise 
de façon essentielle l’interpolation multi-dimensionnelle. 

2. L’interpolation par les fonctions splines : Elle est plus stable que 
l’interpolation polynomiale, est largement utilisée dans tous les pro- 
grammes de dessin assisté par ordinateur, conception assistée par ordi- 
nateur ou plus généralement de graphisme. 

3. Les séries de Fourier et leur analogue discret, la transformation de 
Fourier discrète : Elles sont un moyen très utile pour l’approximation 
des fonctions périodiques. 

Remarque : L’analyse de Fourier est à la base de nombreuses applications, 
par exemple en traitement du signal. 

Remarque : Une façon naturelle d’approcher les fonctions périodiques est 
d’utiliser les polynômes trigonométrique. 

Nous allons nous limiter à l’introduction de l’interpolation Polynomiale : 
c’est la façon la plus classique et la plus simple d’approcher une fonction. Elle 
consiste à déterminer un polynôme P n (x) de degré n qui puisse remplacer 
lors des applications la fonction f(x). 

De plus, c’est un outil efficace pour : 

• Calculer, pour x donné, une approximation de f(x) en calculant P n (x) 

• Construire : 

1. des méthodes d’intégration numérique 

2. des méthodes de différentation 

3. des méthodes d’approximation des équations différentielles 
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4 . ... 


(nous reviendrons en détails sur ces points dans les chapitres suivants). 

Le principe est simple, le procédé est le suivant : 

• On choisit (ou on se donne) (n + 1) points x 0 ,Xi, ...,x n . 

• On calcule y 0 = f{x 0 ), ...,y n = f{x n ) 
ou on se donne (x i7 yi), i — 0 , ..., n . 

• On cherche un polynôme de degré n tel que P n (x t ) = y i, i — 0, ..., n . 

Remarque : 

1) Les points (^,2/;);= o, n sont- appelés points d’interpolation. 

2) Si la fonction / est connue seulement par ses valeurs en quelques points, 
les (n + 1) points xq,xi, ...,x n sont fixés.. 

3) Si on veut que P m (xi) = f(xi) et P' m (xi ) = f(xi), i = 0, ...,n , on 
obtient l’interpolation dite d’Hermite. 

La notion d’interpolation polynomiale est la façon la plus simple d’obtenir 
une telle approximation. 

Nous allons montrer l’existence d’un tel polynôme P n (x) = a n x n + ... + oo 
en le construisant effectivement. 

Il existe plusieurs techniques pour calculer P n (x). Les plus connues sont 
celles de Lagrange et de Newton-Côtes. Elles produisent en fin de compte le 
même résultat. Chaque méthode a ses avantages et ses inconvénients. 

Nous allons en fait le faire des deux façons : 

1. Une méthode directe basée sur la résolution d’un système linéaire 

2. Une méthode itérative due à Lagrange. 

Nous terminerons ce chapitre par : 

1. Une brève discution sur l’erreur d’interpolation polynomiale 

2. Une brève description du principe de la méthode itérée de Newton- 
Côtes 
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4.2 Une méthode directe basée sur la résolu- 
tion d’un système linéaire: 

• On se donne (n + 1) points x 0 , Xi, ...,x n . 

• On calcule y 0 = f(x 0 ), ...,y n = f(x n ) . 

• On cherche un polynôme de degré n tel que P n (xi ) = y,, i — 0, n . 

Écrivons explicitement P n (xi ) = y, . 

a n x™ + a n _i£" _1 + ... + aiXi + a 0 = y i: i = 0, ..., n 

On peut réécrire ces ( n + 1) équations sous forme matricielle : 


( Xq Xq 1 • • • Xq 1 \ 


^ On ^ 



X 1 X l~ l " ' X 1 1 


On— 1 

— 

yi 

ryP rpH 1 ... /y* 1 ; 

\ ^ n ^ n x / 


\ Oo / 


\ y n ) 


La matrice de ce système est une matrice de type Vandermonde. 

On montre que son déterminant est 

det = — Xj ) 

i<j 

On a det ^ 0 si tous les x r sont distincts. On peut donc trouver un unique 
vecteur de coefficients (a n , ...,a o) résolvant le problème. 

Il est connu (à admettre) que les matrices du type Vandermonde devien- 
nent très mal conditionnées lorsque n augmente (elle sont très sensible aux 
erreurs d’arrondies). 

Dans la pratique, cette méthode n’est à utiliser que si n < 3 . Il serait à 
la fois inutile et dangereux de vouloir l’utiliser pour n grand. 
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4.3 Une méthode itérative : Méthode de La- 


grange 

4.3.1 Interpolation Linéaire : 

On considère deux points (x 0 ,y 0 ), (. Xi,yi ) avec : 


f X 0 ^ X X 

\ Vo = f(x o) et y x = f(x i). 


Pour déterminer le polynôme P\{x) de dégré 1 (d’équation : y = ax + b) 
qui passe par deux points distincts (xoj2/o)j {xi,yi) (xo ^ X\). On peut: 

1) Résoudre le système d’équations: 


f ax o + b = y 0 
1 ax\ + b = yi 


d’où 


On a : 


et 


a 

b 


(yi-yo) 


(xi—Xo) 


y o - ax o = 


x iyo •Eoyi 
Xl-XO 


Pi{x) 


(j/i ~ Vo ) x + , xiyp - xpyi 

(xi - X 0 ) X ! - X 0 


Pi(x 0 ) = Vo et Pi(xi) = yi 


2) Poser 


On a: 


L 0 (x) 

Li(x) 


X — Xi 

Xo - Xi 
X — Xq 

X\ - x 0 


Lk (Xi) 


0 si i k 

1 si i — k 
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Ainsi, 


Pi(x) = yoL 0 (x ) + yiLi(x) 

(x-xi) X~Xq 

y ° (x 0 - Xi) 1/1 Xi - x 0 
= 0 j/l - I/o) T + , -Eij/o - Sq2/1 v 
(xi-x 0 ) Xi-Xo 

On a : 


Pi(x 0 ) = I/o et A(xi) = l/l 


car 


Lk (xj) 


0 si i ^ k 

1 si i — k 


Ces deux procédés déterminent évidemment le même polynôme de dégré 
1 (la même droite). 

Si maintenant, on veut déterminer le polynôme de degré 2 qui passe par 
trois (3) points distincts alors: 

i) la première expression de I\ (x) est inadéquate (il faut refaire les 

calculs) 

ii) la dexième expression se prête assez facilement à une généralisa- 
tion par récurrence. 

Exemple : 

Déterminer le polynôme d’interpolation -Pi(x) de degré 1 tel que 
-Pi(xj) = /(x*), i = 0, 1 

avec y t = ffa ) i = 0, 1 , (x 0 , y 0 ) = (0, 1), (xi, yi) = (2, 5) 

On a déterminé le polynôme d’interpolation qui passe par les 2 points : 
(0,1) et (2,5) 

D’après la méthode de Lagrange, 


Pi(x) = y 0 L 0 (x ) + yiLi(x) 

(x - Xi) x- x 0 

= yo— — - + vn- 


(x 0 - Xi) 


Xi - x 0 


^ 1 ( x ~ 2 ) | 5 (^ — 0) 


(0-2) 

= 2x + 1 


(2-0) 
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4.3.2 Interpolation parabolique 

On considère trois points (xo,2/o), (. Xi,yi ) et (x 2 ,y 2 ) avec : 


( x 0 xi ,et x 0 ^ x 2 et x x ^ x 2 
1 yo = f(x 0 ), yi = f(x 1 ) et y 2 = f(x 2 ). 


Pour déterminer le polynôme P 2 (x) de dégré 2, d’équation y = ax 2 +bx+c 
qui passe par trois points distincts (xo, yo), (aq, yi) et (x 2 ,y 2 ), il suffit de 
poser: 


On a : 


Ainsi 


(x — Xi){;x — x 2 ) 
(x 0 - Xi)(x 0 - x 2 ) 
(x - x 0 )(a; - x 2 ) 

Oo - £ O )0 d - x 2 ) 
(x — Xo)(x — X\) 
(x 2 - x 0 )(x 2 - Xi) 


-b'k {Xi) 


0 si i 7 ^ k 

1 si i = k 


Pi{x) = yoL 0 (x ) + 2/1 Z/l (x) + y 2 L 2 (x) 

(x — Xi)(x — x 2 ) (x — x 0 )(x — x 2 ) {x — Xq){;x — X\) 

y °(x 0 - aq)Oo - x 2 ) 2/1 (xi - x 0 )(x 1 - x 2 ) Vl (x 2 - x 0 )(x 2 - x x ) 

est le polynôme d’interpolation polynomiale associé. 

Exemple : 

Déterminer le polynôme d’interpolation P 2 {x) de degré 2 tel que 

Pt.(x%) = f(xi), i = 0,1 et 2 

avec y, = f(x t ) i = 0, 1 et 2 , (x 0 , y 0 ) = (0, 1), (x 1 ,y 1 ) = (1, 2) et (x 2 , y 2 ) = 
(2,5) 

On a déterminé le polynôme d’interpolation qui passe par les 3 points : 
(0,1), (1,2) et (2,5) 

D’après la méthode de Lagrange, 
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p*{x) = yoL 0 (x ) + ViL^x) + y 2 L 2 (x) 

(x — X\)(x — x 2 ) (x — Xo)(x — x 2 ) (x — Xq)(x — X\) 

y °(x 0 - Xi)(x 0 - X 2 ) Vl {x x - x 0 )(xi - x 2 ) V2 {x 2 - x 0 )(x 2 - Xi) 

= A x ~ 1)Qk-2) (x)(x- 2) jx)(x- 1) 

(-l)(-2) (1)(-1) (2)(1) 

= x 2 + 1 

Remarque : 

1) Pour calculer P 2 (x) ,ou n’a pas utilisé le polynôme I\ (x) calculé dans 
l'exemple précédent et pourtant on avait deux points communs. 

2) Li(x),i = 0, 1, 2 sont des polynômes de degré 2 : 

Lo(x) = = \{x - 1) (x - 2) = 2a; 2 - \x + 1 

Li(x) = = —x (x — 2) = — x 2 + 2x 

t (, r \ _ Od(z-i) _ î ™ (, r p — i~2 _ î 

1J 2\ d - ) (2) (1) 2 X y L D 2 J ' 2 X 

On considère (Lj(x)) i= 0)2 comme une base de l’interpolation polynomiale 
quadratique 

Dans l’intervallle [0,2], il existe plusieurs fonctions f{x) qui passent par 
les 3 points (xo,yo) = (0, 1), (aq, yi) = (1,2) et (. x 2 ,y 2 ) = (2,5) mais elle ne 
seront pas approchées par Poix) = x 1 + 1 de la. même façon. 


4.3.3 Interpolation de Lagrange 

• On choisit n + 1 points Xo,Xi, ...,x n . 

• On calcule y 0 = f(x 0 ), ...,y n = f(x n ) . 

• On cherche un polynôme de degré n tel que P n (xi ) = y,, i — 0, ..., n . 


On introduit les coefficients d’interpolation de Lagrange. 


(X ~ X 0 )...(X - gfc-i)(g - x k+1 )...(x - x n ) 
(x k - x 0 )...(x k - æ fe _i)(x fc - x k+1 ...(x k - x n ) 


L k (x) 


n 


3=0 j^k 


(; X — Xj ) 
(X k ~ Xj) 
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L k (x) est un polynôme de degré n, 


P k (.Xi) 


0 si i k 

1 si i = k 


Donc 

n 

p{x ) = y 0 L 0 (x) + yiL^x) + ... + y n L n (x ) = ^ y k L k (x ) 

k=0 

est un polynôme de degré n qui vérifie bien P(xi ) = y* 

Propriété : Le Polynôme d’interpolation polynomiale est unique. 

En effet si P(x) et Q(x) sont deux polynômes d’interpolation alors : 
P(x) — Q(x) est un polynôme de degré n pour lequel 


P(xi) - Q{xi) = 0, i = 0,...,n. 


Ce polynôme de degré < n ayant n + 1 racines, il est identiquement nul. 


Exemple : 

On suppose que f(x) = ffx et que (x 0 ,y 0 ) = (0, 0), (aq, yi) = (1,1) et 

(x 2 ,V2) = ( 8 , 2 ) 

1) Déterminer le polynôme P 2 (a;) d’interpolation polynomiale qui passent 
par les points (x h yi ) i= 0)2 

On a à déterminer le polynôme d’interpolation qui passe par les 3 points 
: (0,0), (1,1) et (8,2) 

D’après la méthode de Lagrange, 


P 2 {x) = y 0 L 0 (x) + yiL^x) + y 2 L 2 {x ) 

(x — xi)(x — x 2 ) (x — x 0 )(a; — x 2 ) (x — Xq)(x — xP 

= yoj- -TT- -T + y 17 - TT7- ITT + 


(x 0 - Xi)(x 0 - X 2 ) (xi - x 0 )(xi - X 2 ) (x 2 — x 0 )(x 2 — Xi) 

(x - l)(x-2) (x-0)(x-8) {x- 0)(a;- 1) 

(0-l)(0-2) (l-0)(l-8) (8 — 0)(8 — 1) 

31 


P 2 (x) = x H X 


28 


28 


On a bien P 2 (0) = 0,P 2 (1) = 1 et P 2 (8) = -J^(8) 2 + f|8 = 2 
2) Calculer P 2 (x) et f(x) = <ÿx pour x = 0.5, 0.95, 1, 1.5 et 3. Conclusion. 
On a : 
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X 

f(x) 

p 2 = ~è x ‘ 2 + ë x 

0.5 

0.793 7 

0.526 79 

0.95 

0.983 05 

0.955 09 

1 

1 

1 

1.5 

1.144 7 

1.4196 

3 

1.442 2 

= 2.3571 


L’interpolation polynomiale de degré 2 ne fournit de résultat acceptable 
qu’au voisinage des points d’interpolation ici 1. 

3) Tracer le graphe de /(x) et P 2 (x) .Conclusion. 

On voit que dans l’intervalle [2,6] , P 2 (x) fournit une mauvaise approxi- 
mation de f(x). 

Pour x donne, P 2 (x) fournira une bonne approximation de /(x) si x est 
voisin de 0, 1 et 8. 

Remarque : 

1) En pratique, on utilise l’interpolation polynomiale avec des polynômes 
de dégré n assez grand ou l’interpolation polynomiale par morceaux. Ainsi 
dans l’exemple précèdent, il faut augmenter le nombre de points d’interpolations. 

2) Si les valeurs ijk sont des valeurs expérimentales. L’interpolation poly- 
nomiale est une technique peu appropriée pour de telles situations. Les 
polynômes de degré élevé sont sensibles à la perturbation des données. 

3) La méthode de Lagrange s’adapte mal au changement du nombre de 
points (x u yi){. On ne peut utiliser les coefficients de Lagrange si on passe 
de n à ( n + 1) points. 

4) Phénomène de RUNGE (fonction de Runge) : L’interpolation 
polynomiale ne fournit pas une bonne approximation de la fonction /(x) = 
1+ J 5j2 . Si on augmente le nombre de points d’interpolation le résultat devient 
plus mauvais. (A admettre). 


4.4 Interpolation Itérée de Newton-Côtes 

• On choisit n + 1 points xo,Xi, ...,x n . 

• On calcule y 0 = /(x 0 ), ...,y n = f(x n ) . 

• On cherche un polynôme de degré n tel que P„(x;) — y i, i — 0, ..., n . 
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L’Interpolation Itérée de Newton-Côtes est un procédé itératif qui permet 
de calculer le polynôme d’interpolation P n (x) de dégré n basé sur (n + 1) 
points (. Xi,yi)i=o, n à partir du polynôme d’interpolation P^ n _i^(x) de dégré 
( n — 1) basé sur n points (xi,yi) i= o,(n-i) , en posant : 

P n {x) = P (n _ i)(x) + C(x), n> 1 

avec 


C(x) 

Q"n 


a n (x - x 0 )(x - xi)...(x - X(„_i)) 


n 


z 


f(Xk) 

( Xk Xq) ...{Xk 1 *^(fc+l) ) • • • 


Xn) 


Les coéfficients a n sont appelés différences divisées d’ordre n de la fonction 
/, on note : 


(Xn f [y'Cb X n ] 

• On appelle "différence divisée d’ordre 0 de / en un point x” la valeur 
définie par 

/ N = f(x) 

• Différence "divisée d’ordre 1 de / en deux points x et y ” la valeur 
définie par 


/ [x, y } 


f[x\-f [y] 

x-y 


on a 


/ [x, y] 


f(x) + f(y ) 

x — y y — x 


• Différence "divisée d’ordre 2 de / en deux points x,y et z” la valeur 
définie par 


/ [x, y, z 


f [x, y}- f [y, z 

f(x) + f(y) + f(z ) 

(■ x -y)(x- z ) (y -x)(y- z ) (z - x) (z - y) 
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et plus généralement: 


n 

f[x 1 ,X 2 ,...,X n ] = ^ 
i= 1 


f(Xi) 


n 

k = 1 
k 7^2 


{Xi 


Xk) 


Remarque: 

Les différences divisées sont indépendants de l’ordre des points. 


Quel est le lien entre f(x) et lex différences divisées? 

Soit x un point autre que les n + 1 points x i: i = 1, n. On a 


mais comme 


alors 


/ [ x , x 0 \ = 


f(x)-f ko] 


X — Xq 


d! où 


f(x ) = / ko] + (x - Xq) f [x, X 0 ] 


/ [x,X 0 ,Xi] = 


f [x,X 0 ] - f ko, Xi] 


X — Xi 


f(x) = f ko] + (x - Xq) f ko, an] - (x- x 0 )(x - xjf [x, X 0 , Xi] 
en continuant ainsi de proche en proche on obtient: 

f(x) = f ko] + (x - X 0 ) / ko, 2u] + ••• + (x ~ X 0 ) ... (x - X n -i) f [x 0 , X n \ + 

(x - æ 0 )...(x - x n )f [x,x 0 , ...,x n \ 

on vérifie que 

f(x) = P n (x) + L(x)f k, X 0 , ...,X n \ 

où P n (x) est un polynôme de degré n tel que P n (xi) = f(xi), pour i = 0, n. 

C’est donc le polynôme d’interpolation de Lagrange, on l’appelle le polynôme 
de Newton. 

Comme signalé dans l’introduction, l’interpolation polynomiale sera util- 
isé comme outil d’approximation (pour la construction des méthodes d’intégration 
numérique ou des méthodes de dérivation numérique ou des méthodes d’approximation 
des équations différentielles), il est donc fondamental de connaître une ex- 
pression de l’erreur d’interpolation. 
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4.5 Erreur d’interpolation polynomiale : 

L’erreur commise lors d’une interpolation est une question fondamentale en 
analyse numérique: 

• elle renseigne à priori sur la nature de cette erreur 

• elle fournit des informations sur les termes qui y participent 

• elle permet d’avoir un ordre de grandeur de l’erreur commise. 

Nous allons énoncer un résultat qui répond à ces interrogations dans le 
cas où la fonction / est régulière (de classe C p ,p assez grand). 

Théorème : 

Soient / une fonction de classe C n+l dans / et , (xj) i= 0 , n (n + 1) points 
distincts dans / avec xq < x± < ... < x n 

Alors pour tout x G [xo,Xn\ , il existe ( = ((x) tel que 


f( n + i)(7ù fO+Of/q 

f(x ) - Pn(x) = - x 0 )(x - X!)...(x - X n ) = ——^ L (x) 

où 



P n {x 

avec 

L k (x 

et 

L{x 


n 

yoLo(x ) + yiL^x) + ... + y n L n (x ) = ^ ykL k (x ) 

k = 0 


n 


n 

3 = 0 j+k 


(x - Xj) 
(■ Xk - Xj) 


(x - x 0 )(x - x 1 )...(x - x n ) 


Pjx) est le polynôme d’interpolation de Lagrange. 

Remarque : 

1) Cette formule montre que : 

i) l’erreur est nulle pour x = Xi i.e. x est un point d’interpolation. 

ii) l’erreur dépend de la fonction considérée ( de f (n+1 >) et des points 
d’interpolations (xj)j. 

2) Cette formule d’erreur permet de trouver des formules d’erreur pour 
l’intégration numérique et la différentiabilité numérique. 
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Dans le cas de l’erreur d’interpolation à partir de la forme de Newton, on 
a: 

f(x) - P n (x) = L{x).f[x,x o, x n \. 

Comme on a la même fonction / selon les mêmes points x, pour i = 0, n. 
il s’agit de deux formes du même polynôme, et l’erreur d’interpolation est la 
même, d’où 


f(x) - P n (x) 


f( n + 1)(0 

( n + X )l L ( x ) = L (x).f[x,x 0 ,...,x n ]. 


Exemple : 

Déterminer l’erreur d’interpolation polynomiale dans le cas de l’interpolation 
parabolique 

On approche la fonction f(x) par la parabole passant par les points 
(zo,î/o) = (0, 1), (a?i, 2 /i) = (1,2) et (x 2 ,y 2 ) = (2,5). 

Le polynôme d’interpolation P 2 (x) de degré 2 tel que P 2 (xi) = /(ay), 
i = 0, 1 et 2 

avec t/j = /( x t ) i = 0, 1 et 2 , (x 0 , y 0 ) = (0, 1 ), , y, j = (1,2) et (x 2 , y 2 ) = 

(2,5) 

D’après la méthode de Lagrange, 


Pi{x) = y 0 L 0 (x) + yiLi(x) + y 2 L 2 (x) 

= (x-l)(x-2) (x)(x-2) (x)(x-l) 

( 1 ) ( 2 ) + ( 1 )(- 1 ) + ( 2 )( 1 ) 

= x 2 + 1 

D’après le théorème précédent, 


f(x) - P 2 (x) 


/( 3)(0 

— — — ( X - x 0 ){x - x x )(x - x 2 ) 

f m (Q , ,w 

~——x{x - l)(x- 2) 


Si f'( 3 \x) < M alors 
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Oi| h-* 


Va; G [0,2] , \f(x) - P 2 (x)\ 

M 

I f(x) - Pi{x)\ < — \x(x - l)(x - 2)| 
o 

M 

< — x(x — l)(a; — 2) 

6 

< 6.4 * 10^ 2 * M. 

(le maximum de u(x) = x(x — l)(x — 2) est atteint en x* = ; d’où 

u(x*) = - l) - 2 ) = 0.0 641 5 ~ 6.4 * 10” 2 ). 

4.6 Exercices: 
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Série Interpolation Numérique 


Exercice I : 

1) Déterminer par une méthode directe basée sur la résolution d’un sys- 
tème linéaire, le polynôme d’interpolation P\ (x) de degré 1 tel que Pi(xj ) = 
/(x*), i = 0, 1 avec y { = f(xi ) i = 0,1, (x 0 , Vo) = (-2, 4) et (xi, yi) = (2, 8) 

2) Déterminer par une méthode directe basée sur la résolution d’un sys- 
tème linéaire, le polynôme d’interpolation P 2 (x) de degré 2 tel que P 2 (xj) = 

i = 0, 1 et 2 avec y, = f(x t ) i = 0, 1 et 2, (x 0 ,y 0 ) = (-2,4), (xi,yi) = 
(0,2) et (x 2 ,y 2 ) = (2,8). Conclusion. 

Exercice II : 

1) Déterminer par la méthode de Lagrange, le polynôme d’interpolation 
Pi(x) de degré 1 tel que Pi(xj) = /(x*), i = 0, 1 où y t = /(x* ) i = 0,1, 
(zo,î/o) = (-2,4) et (xi,yi) = (2,8) 

2) Déterminer par la méthode de Lagrange, le polynôme d’interpolation 
P 2 (x) de degré 2 tel que P 2 (x;) = /(xj), i = 0, 1 et 2 où yi — /(x; ) i = 0, 1 
et 2, (x 0 ,î/o) = (-2,4), (xi,t/i) = (0,2) et (x 2 ,y 2 ) = (2,8) 


Exercice III : 

1) Déterminer par la méthode de Newton-Côtes, le polynôme d’interpolation 
Pi(x) de degré 1 tel que Pi(xj) = /(x*), i = 0, 1 où y % = /(x*) i = 0,1, 
(®o,yo) = (-2,4) et (xi,j/i) = (2,8). 

2) Déterminer par la méthode de de Newton, le polynôme d’interpolation 
P 2 (x) de degré 2 tel que P 2 (x*) = /(x*), i — 0, 1 et 2 où yi — f{xi ) i = 0, 1 
et 2, (x 0 ,yo) = (-2,4), (xi,yi) = (0,2) et (x 2 ,y 2 ) = (2,8). Conclusion. 


Exercice IV : 

On suppose que (x 0 , y 0 ) = (0, 0), (x x , j/i) = (1, 1) et (x 2 , y 2 ) = (2, 8) 

1) Déterminer par la méthode de Lagrange, le polynôme d’interpolation 
P 2 (x) de degré 2 tel que P 2 (x*) = y*; i = 0, 1, 2. 

2) Tracer le graphe des fonctions P 2 (x) = 3x 2 — 2x et /(x) = x 3 dans 
l’intervalle [0,2] . 

3) Calculer P 2 (x) et /(x) = x 3 pour x = 0.9, 1.1, 1.99, 2.1 et 5. Conclusion. 

4) Déterminer l’erreur commise si on en remplace dans l’intervalle [0,2] , 
/(x) = x 3 par P 2 (x) = 3x 2 — 2x. 

Exercice V : 

On suppose que (x 0 ,y 0 ) = (0, 1), (xi, yi) = (0.5, es), (x 2 ,y 2 ) = (l,e) 
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1) Déterminer par la méthode de Lagrange, le polynôme d’interpolation 
Pïix) de degré 2 tel que P-z(xi) — y^i — 0, 1, 2 et 3. 

2) i) Déterminer une expression de l’erreur d’interpolation polynomiale, 
ii) Déterminer une borne de l’erreur d’interpolation polynomiale. In- 
dépendantes de £ où £ = £(x). 

ii) Déterminer une borne de l’erreur d’interpolation polynomiale. In- 
dépendantes de £ et de x. 

i 


X S. El Bernoussi, S. El Hajji et A. Sayah 


54 



Chapitre 5 

Intégration et dérivation 
numérique. 


5.1 Introduction : 


Si / est une fonction dérivable sur [a, b], la dérivée en c G ]a,b[ est définie 
par: 


f(c) 
où A /(c) 


linr 

h-> 0 


A /(c) 

h 


f(c + h) - /(c) 


Si / est une fonction continue sur [a, b\, l’integrale de / sur [a, b] est définie 

par 


f(x)dx 


fini R(h) 

h^O 


où R(h) = f(a + kh).h 

k = 1 


R(h) est la somme de Riemann avec h = 

On sait déterminer f { c) ” exactement” pour / définie à partir de fonctions 
élémentaires (exp: sinx, e x , ln.x, ...). 

On sait aussi calculer f'’ f{x)dx en utilisant les théorèmes fondamentaux 
d’intégration pour une fonction continue sur [a, b ], et on a f b f(x)dx = F (b) — 
F (a) où F(x) est une primitive de f(x). 
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Mais il existe des fonctions très simples comme ou \/ cos 2 x + 3 sin 2 x 

qui n’ont pas de primitive connue, donc, comment peut-on intégrer de telles 
fonctions entre a et b? 

D’autre part / peut-être connue seulement en quelques points et sa for- 
mule est inconnue (exp: résultats experimentaux,...), donc comment peut-on 
dériver ou intégrer ses fonctions? 

Du point de vue numérique, la solution à ce problème est immédiate: nous 
avons vu, dans les chapitres précédents, comment approximer une fonction 
par une fonction plus simple, facile à dériver ou à intégrer. 

De façon précise si P(x) est une approximation de / dans l’intervalle [a, b], 
nous nous proposons d’étudier les approximations: 

f(y) ~ P' (y) y e k b] 

et 

pb pb 

/ f(x)dx ^ / P(x)dx. 

J a J a 


5.2 Dérivation. 

La dérivation numérique nous permet de trouver une estimation de la dérivée 
ou de la pente d’une fonction, en utilisant seulement un ensemble discret de 
points. 


5.2.1 Dérivée première. 

Soit / une fonction connue seulement par sa valeur en (n + 1) points donnés 
Xi i — 0, 1, ..., n distincts. 

Les formules de différence les plus simples basées sur l’utilisation de la 
ligne droite pour interpoler les données ulilisent deux points pour estimer la 
dérivée. 

On suppose connue la valeur de la fonction en Xi-i,Xi et Xi + q on pose 
f(xi- 1 ) = Vi- 1 , f{xi) = y t et f(x i+ 1 ) = y i+ 1 . 

Si on suppose que l’espace entre deux points successifs est constant, donc 
on pose h = Xi — aq_i = x i+ i — x, L . 

Alors les formules standarts en deux points sont: 

Formule de différence progressive : 
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f(Xi+l) - f(Xi) 


%i -\- 1 %i 

Formule de différence régressive : 


f(Xi) 


f(Xi) ~ f{Xi- 1) 


%i — 1 

Formule de différence centrale 


V %+ 1 - Vi 


%i + 1 


Vi - Vi - 1 
1 


w/ x /(®i+i) - /te-i) y*+i - î/*-i 

J (Xi) « = . 

•^i+1 *£i— 1 -^i+1 Xf— i 

Les trois formules classiques de différences sont visualisées sur la figure 
suivante, et sont les conséquences de la définition de la dérivée: 

Exemple : 

Pour illustrer les trois formules, considérons les données suivantes: 
(xo,Vo) = (1,2);(xi,î/i) = (2,4 )\(x 2 ,y 2 ) = (3, 8); (x 3 ,y 3 ) = (4,16) et 
(x 4 ,y 4 ) = (5,32). 

Nous voulons estimer la valeur de f'(x 2 ). 

Progressive: f'(x) ~ = 8. 

Régressive : f'(x) ~ = fzl = 

Centrale : /'(x) « = 6. 

Les données ont été calculé pour la fonction /(x) = 2 X . f(x) = 2 X ln(2) 
et pour x = 3 /'( 3) = 2 3 ln(2) = 5.544. 


Remarque: 

Les formules de différences classiques peuvent être trouvées en utilisant 
la formule de Taylor. 

f(x + h) = /(x) + hf(x) + y/" (y), 
x < rj < x + h 

• Formule progressive: 

h Xî -\- 1 x^ 

/'(,.) = /fe) - in,) 

Xi<rj < x i+1 


57 



l’erreur est |/ ,/ ( 7 ?) donc en O (h). Cette formule peut être trouvée aussi 
en utilisant le polynôme d’interpolation de Lagrange pour les points 

(xi,f(xi)) et (x i+1 , f(x i+1 )). 

• Formule régressive: 


h - Xi - i 

/'M = M z^A + 'i rM 

Xi - 1 < rj <Xj 


La formule de différence centrale de la dérivée en x t peut être trouvée en 
utilisant la formule de Taylor d’ordre 3 avec h = Xj+i — Xi = Xi — Xi-\ 

f{Xi+ 1) = f(Xi) + hf'{Xi ) + y/"(Xj) + i) 

f{xi- 1 ) = f{xi) - hf(xi) + y/"(xi) - ^f"\r] 2 ) 

Xi < r/i < x i+ i, Xi_i < r/2 < x. 


si on suppose que f" est continue sur [x*_ i,x i+ i] on peut écrire la formule 
suivante: 


< V < Xi+i 


l’erreur est donc en 0[h 2 ). La formule de différence centrale peut 

aussi être trouvée à partir du polynôme d’intérpolation de Lagrange en 3 
points. 


On peut interpoler les données par un polynôme au lieu d’utiliser la 
droite, nous obtenons alors les formules de différence qui utilisent plus de 
deux points. On suppose que le pas h est constant. 

Formule de différence progressive utilisant trois points: 

~/(xj +2 ) + 4/(xj+i) - 3 f{xj) 

Xi- 1-2 Xi 

Formule de différence régressive utilisant trois points: 


f\Xi) 
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f{Xi) 


3 f(xi) - 4/(xi_i) + f(xi- 2 ) 


%i— 2 

Exemple : Formules de différence en trois points: 


En utilisant les données de l’exemple précédent, on trouve: 

-PU \ ^ — 32 + 4 ( 16 )— 3 ( 8 ) A 

j [Xi) ~ — — — 4 progressive. 

f(xi ) « 3 ( 8 )~F'P +2 _ g régressive. 


5.2.2 Formule générale en trois points. 

La formule d’approximation en 3 points de la dérivée première, basée sur le 
polynôme d’interpolation de Lagrange, n’utilise pas des points équidistants. 

Etant donné trois points (xi,yi)] (£2,1/2) e t (£3,1/3) avec £1 < £2 < £3, la 
formule suivante permet d’approcher la dérivée en un point x G [£+, £3] . Les 
dérivées aux points x t sont les suivantes: 


f'M 

f'M 

f'M 


2 x 1 - £2 - £3 , £1 - £3 , £1 - £2 

(£1 - £2) (£1 - £3 ) Vl (£2 - £1) (£2 - £3 ) Vl (£3 - £1) (£3 - £2) 2/31 
£2 - £3 , 2£ 2 - £1 - £3 , £2 - £1 

(£1 - £2) (£1 - £3 ) V1 (£2 - £1) (£2 - £3 ) V2 (£3 - £1) (£3 - £2) 2/31 
£3 - £2 , £3-^1 , 2£ 3 - £2 - £1 

(£1 - £2) (£1 - £3 ) V1 (£2 - £1) (£2 - £3 ) V2 (£3 - £1) (£3 - £2) 2/31 


Le polynôme de Lagrange est donnée par 


P(x) = L 1 (x)y 1 + L 2 (x)y 2 + L 3 (x)y 3 
où 

(x - £ 2 )(£ - £ 3 ) 

1 X (£1 - £ 2 )(£i - £3) 

J , A = (£ - Xi)(x - £3) 

2 ' 1 (x 2 — X]_)(x 2 — £ 3 ) 

(x - £i)(£ - £2) 

3 X (£3 - £1) (£3 - x 2 ) 

L’approximation de la dérivée première est donnée par f'(x) ~ P'(x), qui 
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peut s’ecrire 


donc 

/'(*) 


P'(x ) = L , 1 (x)y 1 + L' 2 (x)y 2 + L' 3 (x)y 3 
où 


L[(x) 
L' 2 (x ) 
L 3 { x ) 


2x — x 2 — x 3 
{xi - x 2 )(xi - x 3 ) 
2a; — x\ — x 3 
(x 2 - Xi){x 2 - x 3 ) 
2a; — x\ — x 2 
(x 3 - xi)(x 3 - a; 2 ) 


2a; — x 2 — x 3 2x — x\— x 3 

+ 77 17777 771/2 + 


2a; — x\ — x 2 


(xi - x 2 )(xi - x 3 ) (x 2 — x\){x 2 — x 3 ) (x 3 — X]_)(x 3 — x 2 ) 


2 / 3 - 


5.2.3 Dérivées d’ordre supérieur. 

Les formules de dérivées d’ordre supérieur, peuvent être trouvées à partir des 
dérivées du polynôme de Lagrange ou en utilisant les formules de Taylor. 

Par exemple, étant donné 3 points x i: x i+ i équidistants, la formule 
de la dérivée seconde est donnée par: 

f"( x i ) = jp[f(Xi+ O - 2 f(Xi) + f(Xi- 1)] 
l’erreur est en 0[h 2 ). 

Dérivée seconde à partir du polynôme de Taylor. 
f(x + h) = /( x) + hf(x) + Ç/"(x) + ^f"(x) + ÿ/ (4) (ry i) 
f(x- h) = f(x)~ hf(x ) + Ç/"(x) ~ f /"'(*) + ÿ/ (4) 0? 2 ) 
x < rji < x + h et x — h < r /2 < x. 

f /// n _ f( x + h) + f(x -h)- 2 f(x) 

J W - h 2 

l’erreur est en 0[h 2 ). 

Pour obtenir les formules de la troisième et la quatrième dérivée, on prend 
une combinaison linéaire des développement de Taylor, pour f(x + 2h), f(x + 
h), f(x — h) et f(x — 2 h). 
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La table suivante donne différentes formules centrales toutes en O (h 2 ): 


ffa) - ^ [f(Xi+ 1) - f(Xi- 1 )] 
f"(Xi) ~ 1 ) “ 2 f( x i) + f( x i~ l)] 

f"(xi) ~ 2^3 [/( x *+2) - 2 /(x i+ i) + 2/ (xj_i) - / (a?i_ 2 )] 

/ (4) (^i) ^ ^ [/(x» +2 ) - 4/(x i+ i) + 6/(xi) - 4/(xj_i) + f(xi - 2 )] • 

En utilisant les polynômes d’interpolation de Lagrange les dérivées d’ordre 
p sont calculées par: 


f {p \a) ~ J ZAi(ai)f(xi ) 

i = 0 

OÙ 

Ai (a) = L[ p \a) p <n 

n 

Aj(o:)x* = 0 0 < A; < p — 1 

i=0 

n 

Ai(a)Xi = k(k — 1 )...(k — p + l)a fc_p p < k < n. 
i = 0 

Remarque : 

La fomule est exacte pour les polynômes de degrés < n. 

Le système linéaire donnant les A^a) a un déterminant de type Vander- 
monde différent de zéro si les Xi sont distincts. 

Les A t (a) sont indépendants de / et peuvent être calculés une fois pour 
toutes. 
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5.2.4 Etude de l’erreur commise. 


D’après le chapitre précédent, si / est connue en (n + 1) points Xi, i = 0 , n 
alors f(x) = P n (x ) + e(x), où e(x) est l’erreur d’interpolation. En dérivant 
on obtient 


î'{x ) = P' n (x) + e\x) 


i / \ oe 

et e (x) = — 


X] Ai(x).f(xi ) + e' 

i=0 

d 


X) 


dx \{n + 1)! 

1 


r(z)-/ < ” +1> (îjf = d-{L{x).f[x 0 , -,x„,x\) 


(n + 1) 






On remarque tout de suite que l’erreur de dérivation est nulle si / est un 
polynôme de degré inférieur ou égale à n. Si on prend pour x un point 
Xi, le second terme de la dérnière somme s’annule, sinon il faut connaître 
£ (/ (n+1) (&)) , ce qui est difficile car la fonction x — > £ r étant inconnue. On 
peut donner une forme si / est n + 2 fois dérivable en utilisant la notion de 
différence. En effet 


d 

dx 


(/ (n+1) ( 6 )) 


— (f[x o,...,x n ,x}) 

li /[^Q; .■■,X n ,X + h ] - f[x 0 , ...,X n ,x) 
h^O h, 

lim/[x 0 , ...,x n ,x,x + h] 


lim - 

h-> 01 


f {n+2 \0 x ,h). 


n 


2 )! 


On constate qu’on devra se contenter d’une estimation 


e (x ) |< 


1 

(n + 1)! 


L\x) | M n+ i + 


1 

(■ n + 2 )! 


L(x) | M n+ 2 . 


5.3 Méthodes numériques d’intégration. 

Le but de cette leçon est de calculer numériquement des intégrales définies 
ou indéfinies. Soit / : [a, 6] — > M, une fonction continue donnée. On désire 
approcher numériquement la quantité J ^ f(x)dx. 
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5.3.1 Formules fermées. 

On appelle ainsi les formules quand la fonction f est continue sur lintervalle 
[a, b]. Les points d’interpolation Xi vérifient a = xq < x\ < ... < x n -\ < 
X n = b. 

Formule des rectangles. 

La formule des rectangles est une formule dite à un point xq = a. Le 
polynôme d’interpolation associé est Pq{x) = /(a) et L(x) = x — a pour tout 
x appartenant à [a, b]. D’où 


/(/)~I(P„) = /(<*)(& -a). 

L’interprétation graphique consiste donc à remplacer f b f(x)dx par l’aire 
du rectangle de base [a, b] et de hauteur f(a). 

Formule des trapèzes. 

La formule des trapèzes est une formule à 2 points : xq = a et x\ = b. 
Le polynôme de Lagrange associé à ces deux points est P\{x) = /(a) (yrf) + 

/w(k)-d’où 


/(/) ~ J(Pi) = fpy(x)dx = n " ] ± m {b - a). 

Formule de Simpson. 

La formule de Simpson est une formule à trois points x 0 = a , x i = ^ 
et x 2 = b: . Le polynôme associé à ces trois points est P 2 (æ) = f(a)L 0 (x) + 
f( g ^)Li(x) + f(b)L 3 (x). Notons que 


L 0 (x) 


L i(x) 


(x — X\){;x — b) 
(a — x \ )(a — b) 

(. x — a)(x — b) 
(aq — a)(x i — b) 


L 2 (x) = 


(. x — xi)(x — a) 


( b — xi )(b — a) 
On tire donc la formule suivante: 


L 0 (x)dx 


Li(x)dx 


L 2 (x)dx 


(i b — a) 


4(6 — a) 
6 

(b — a) 
6 ’ 


m * m) 


(i b — a) 
6 


/(«) +4/(^4 ) +m 
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Formules ouvertes. 

On appelle ainsi les formules quand la fonction / est continue sur l’intervalle 
]a,b[. Les points d’interpolation aq vérifient a < Xq < X\ < ... < x n -\ < 
x n < b. 

Formule de Steffensen. 

Il en existe une infinité. 

• Une à 1 points avec Xq = qui donne la formule du milieu suivant: 

/(/) ~ (6- a)/(^y^) 

Cette formule est exacte pour tout polynôme de degré 1. 

• Une à 2 points avec x 0 = et aq = q U j donne la formule 
suivante : 

Cette formule est exacte pour tout polynôme de degré 1. 

• Une à 3 points avec xo = et aq = ^ et aq = qui donne la 
formule suivante : 

m - ^ (v (^) + 2/ (^) - 2/ (^)) . 

Cette formule est exacte pour tous les polynômes de degré 2. 


5.3.2 Etude générale de l’erreur commise. 

Pour que les formules d’intégration numérique données précédement soient 
intéressantes, il faut que l’on puisse estimer l’erreur E(f) = /(/) — I(P n ) 
avec précision. Or si / est suffisamment dérivable, on a 


E(f) = I(f-P n ) 



1 

(n + 1)! 




dx. 


Théorème : Supposons que E(f) = 0 pour les polynômes de degré au plus n 
et que la fonction / G C n+l ([a, 6]). On dit alors que la méthode est d’ordre 
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n+1. Si on pose M n+1 = max | /( n+1 )(x) |, Une première estimation de 

x£[a,b] 

l’erreur est 


E(f) l< 


:M n +i / I L(x) j dx. 


(n + 1 )! 

Théorème : En plus des hypothèses du Th précédent, on suppose que le 
polynôme L(x) ne change pas de signe sur [a, b], alors en utilisant le Th de 
la moyènne pour E(f), on obtient 


E(f) = 


n + 1)! 
V e [a, b ] 


/( n+1 )( rj) / L(x)dx. 


En utilisant ce dernier Théorème on peut estimer les erreurs des méthodes 
vues ci-dessus. 


• Pour la formule du rectangle on a: 

EU) = f'(v) f a ^ x ~ a ) dx = f'iU — 2^ V e [a, b] 

cette méthode est d’ordre 1. 

• Pour la formule du trapèze on a: 

E U) = \ f"U) J a ^ x ~ a ^ x ~ b î dx = ~ 

la méthode de Trapèze est dordre 2. 

• Pour la formule de Simpson on a: 

b — a 5 

2 J ’ 

la méthode de Simpson est dordre 4. 


E(f) = 


fU_M 

90 


Exemple : 

I = e ~ x2 dx, a = 0,Ur 
.36788. 


I ,b = l,/(0) = 1 ,/(§) = .7788, /(l) = 
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1. Rectangle: / ~ /(O) = 1. 

2. Trapèze: I ~ — .68393. 

3. Simpson: I [/( 0) + 4/(4) + /( 1)] = .74718. 

4. La valeur de / à 5 décimales est .74718. 


5.3.3 Formules composées. 

Plutôt que d’augmenter le degré du polynôme d’interpolation, on peut obtenir 
une formule d’intégration en découpant l’intervalle d’intégration en sous- 
intervalles et en appliquant des formules simples sur chacun des sous-intervalles. 
Formule de trapèze. 

Si n est entier, posons 

h= -,Xk = a + kh, k — 0,...,n. 

n 


rb n 1 / p 

I U) = / f(x)dx = ^2( 

J a k = 0 x k 


f(x)dx 


k = 0 L V / 

où rjk e [x k ,Xk+i\ , k = 0, ..., n - 1 

Développant et regroupant les termes qui apparaissent 2 fois, on obtient 


I (f) = ï /(«) + 2^/(a + kh) + f{b) ~ 


En appliquant le Th des valeurs intermédiaires, on peut réécrire l’erreur 
sous la forme 

E(J) = - n 4rw = 
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Ceci nous donne la formule du trapèze composée pour laquelle l’approximation 
est donnée par: 


Tn(J) = \ 


71—1 


/(a) + 2 f(a + kh ) + / ( b ) 


k = 1 


et l’erreur par 

ET(f) = -^d/"(>;)fc 2 . 

Formule de Simpson composée. 

Supposons maintenant que n soit pair, groupant les intervalles 2 à 2 et 
appliquant la formule de Simpson sur x t , a;^], on obtient 


Hf) 


h 

3 


/ (°) +4 ^ 


A: impair 


f(a + kh) + 2 ^ /(a + kh) + f(b) 

k pair 


n f ( 4 ) (v) h s 
2 90 


Ceci nous conduit à la formule de Simpson composée pour laquelle l’approximation 
est donnée par 


-»„(/) 


h 

3 


/(a) +4 ^ 


k impair 


f(a + kh) + 2 ^ /(a + A:/i) + /(6) 

Ai pair 


et l’erreur par 

ES(f) = 

Exemple : Déterminer J^e~ x2 dx. 

Si n désigne le nombre des intervalles utilisés. 


n 

Uf) 

ET(f) 

2 

.73137 

.015 

4 

.74298 

3.84 x 10“ 3 

8 

.74658 

9.58 x 10 -4 

16 

.74676 

1.39 x 10 -4 

32 

.74680 

5.98 x 10“ 5 


Si nous désirons obtenir 6 décimales exactes, il nous faut déterminer h tel 


que 


max | f"{rj) 

y)<7}<\ 


12 


< 5 x 10“ 7 , 


(5.1) 
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Pour une partition régulière xj, = kh,h = donc nous cherchons n tel 


que 


n 2 > — m ax | f'(rj) 


120<r?<l 

or f"(x) = e~ x2 (4x 2 — 2) et f"(x) = e~ x ^4x(3 — 2x 2 ). Puisque f'"(x) ne 
change pas de signe sur [0; 1] , 


5 x 10“ 7 
- 2 ' - „ 2 \ 


max |/"(r / )|=max{|r(0)|,|r(l)|} = 2. 

0 <? 7<1 


On voit que (5.1) sera satisfaite si 

2> io 6 

n >— . 


n > 578. 


Remarque Dans le choix de la précision demandée, il faut tenir compte des 
erreurs d’arrondi et de l’accumutation des erreurs 


5.4 Exercices 
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Série intégration et dérivation. 


Pour les problèmes des exercices (5.4) et (5.4), donner des approximations 
des dérivées dans les cas suivants: 

En utilisant la formule de différence progressive. 

En utilisant la formule de différence régressive. 

En utilisant la formule de différence centrale. 

Exercice : 1 

Approcher y 7 (1.0) si 

a; =[0.8 0.9 1.0 1.1 1.2] 
y = [0.992 0.999 1.000 1.001 1.008] 

Exercise: 2 

1. Approcher y' (4) si 

x = [0 1 4 9 16] 
y=[ 0 1 2 3 4] 

2. Donner une expression de l’erreur de dérivation en x — 4. 

3. Donner une majoration de l’erreur independament de x et de 

Exercise : 3 

Calculer y" (2) si 

x = [0 1 2 3 4] 
y = [0 1 4 9 16] 

Exercise : 4 

Calculer f () 2 sin 2 xdx en utilisant la formule du trapèze et la formule de 
Simpson. Comparer avec le résultat exact. 

Exercise : 5 

Pour le problème PI approcher 1’ intégrale: 
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1. En utilisant la formule de trapèze composée avec 2 intervalles. 

2. En utilisant la formule de trapèze composée avec 10 intervalles. 

3. En utilisant la formule de Simpson avec 2 intervalles. 

4. En utilisant la formule de Simpson composée avec 10 intervalles. 


PI 


x sm(irx)dx 


Exercise : 6 

En utilisant les formules d’estimation d’erreur, trouver les bornes d’erreur 
pour le problème PI dans les cas 1-4, puis calculer la valeur exacte de 
l’integrale et comparer les erreurs exactes ” £’(/)” et les bornes d’erreurs 
trouvées. 

î 


X S. El Bernoussi, S. El Hajji et A. Sayah 
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